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I. INTRODUCTION 


Space radiation poses a serious risk to the health of 
astronauts in space explorations [1]. Space radiation par- 
ticles mainly consist of the galactic cosmic rays (GCR) 
and solar energetic particles (SEP). SEPs can have a very 
large flux in solar particle events such as coronal mass 
ejections but their kinetic energies are typically below a 
few hundred MeV. GCR particles are fully-ionized nuclei 
at much higher average energies, thus they can penetrate 
thick materials and pose serious health hazards, espe- 
cially for long-term human space explorations. Physi- 
cal and biological effects of space radiation are currently 
being studied [2] in order to provide astronauts enough 
protection against space radiation. 

Reliable space radiation transport codes are essen- 
tial for the prediction of radiation environments behind 
shielding and effective reductions of the radiation risk. 
A key input to these transport codes is the physics of 
the interactions between radiation particles and different 
materials. Atomic interactions which cause the ioniza- 
tion energy loss are well known. However, significant 
uncertainties exist in our understanding of the nuclear 
fragmentation processes of different ions. To improve 
the accuracy of the space radiation transport, labora- 
tory experiments have been performed to measure frag- 
ments from the fragmentation of different elemental pro- 
jectiles, e.g., at the NASA Space Radiation Laboratory 
at the Brookhaven National Laboratory. Because space 
radiation particles involve many different ions and they 
fragment inside materials, an important question is what 
partial cross sections affect the space radiation risk more 
and thus need to be better studied. 

Here we present a semi-analytical method in order to 
understand the qualitative features of the effects of nu- 
clear fragmentations on space radiation risks from GCR. 
In Sec.II we describe the semi-analytical method and the 
theoretical results obtained in the thin-shielding limit, 
including a list of the top 30 most important partial 
cross sections for both water and aluminum shielding 
for the 1977 solar minimum GCR environment. Sec.III 
shows our results for thick shielding. In Sec. IV we dis- 
cuss the important effect of the unitarity relation due to 
the baryon number conservation, and provide two simple 
examples on alternate implementation of the unitarity re- 


lation. The section also shows the difference in the sensi- 
tivities between solar maximum and minimum GCR en- 
vironments. We summarize and present our conclusions 
in Sec.V. 


II. THIN-SHIELDING RESULTS FROM THE 
SEMI-ANALYTICAL METHOD 


In the straight-ahead approximation, where fragmen- 
tation products are assumed to keep the same velocity as 
the projectile, the propagation equation is given by [3] 

dJ k {E,x) Jk{E,x) s -' J j(E> x ) 9[uj^E)J k (E,x)} 
dx A k (E) ^Akj(E) + dE 1 


where J k (E,x) represents the flux of particle species k 
with energy E per nucleon at depth x in a material, and 
uj k (E) represents the ionization energy loss. The term 


A k (E) = 


1 

na k (E) 


( 2 ) 


represents the mean-free-path of particle species k in the 
shielding material with n being the number density of 
nuclei in the material and a k {E) the total inelastic cross 
section of particle species k. In the term 


A kj {E) = 


1 

na k j (E) 


( 3 ) 


a k j{E) represents the partial cross section for projectile 
j to produce fragment k. In the thin-shielding limit, i.e. , 
x — > 0, we obtain the following results to first order in x: 


Jk{E,x) 


Jk{E, 0) 


1 + ui' k (E)x + 


J' k (E, 0 ) 
J k (E, 0 ) 


uj k (E)x 


x 

Ak(E) 
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ME, 0 ) 
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( 4 ) 


where ui' k (E) and J' k (E, 0) represent the first derivative 
with respect to energy E. Only the last two terms on the 
right-hand side of Eq. (4) involve the fragmentation cross 
sections and they do not involve the ionization energy loss 
t Jk{E). Therefore, to first order in x, the energy loss in 
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materials does not affect the sensitivity of the shielded 
radiation fluxes to fragmentation cross sections. 

Dose equivalent H is used in this study to relate the 
energy spectra of radiation particles to the radiation risk: 

H = J Q(L)D L dL, (5) 

where Di represents the distribution of absorbed dose as 
a function of L, the unrestricted Linear-Energy- Transfer 
(LET). The International Council on Radiation Protec- 
tion (ICRP) has defined a quality factor, Q{L), which 
roughly summarizes the data on Relative Biological Ef- 
fectiveness (RBE) of different types of radiation relative 
to X-ray or 7-ray radiation. Because RBE depends on the 
biological end-points, and certain types of radiation may 
have qualitatively different effects than X-ray or 7-ray 
radiation, there are significant uncertainties associated 
with the quality factor. In this study we use the quality 
factor defined in the most recent ICRP60 [4]. The dose 
equivalent inside soft tissue, which is simulated by water, 
behind depth x of a shielding material is then given by 

H(x) = —Y,f Jk[E,-x)L h (E)Q(L h (E))dE , (6) 

PT k J 

where pr represents the density of water and L k ( E ) rep- 
resents the LET value in water for particle k at energy 
E per nucleon. 

Eqs. (4) and (6) lead to * 


5H{x) = ^E/4 -LML 3 )Sc j+ ^L k Q(L k )Sa kj 


dE,( 7) 


which relates the uncertainty in the dose equivalent to the 
uncertainty in each partial fragmentation cross section. 

The baryon number conservation relates the total in- 
elastic cross section to the partial fragmentation cross 
sections. When the energy of the projectile is below 
the threshold for baryon-antibaryon productions, about 
6 GeV per nucleon in the laboratory frame, the baryon 
number conservation is expressed as the following unitar- 
ity relation: 

AjffjiE) = (8) 

k 

assuming the straight-ahead approximation, where frag- 
ments of the projectile can be separated from fragments 
of the target material. The uncertainties in the total in- 
elastic cross section and the partial fragmentation cross 
sections are then related via the unitarity relation as 

A j S* j (E) = ^2 A k 6<T k j( E ). (9) 

k 

Therefore not all uncertainties in Eq. (7) are indepen- 
dent. Because GCR particles with kinetic energies from 
about 0.2 GeV to a few GeV per nucleon, i.e., mostly 



FIG. 1: (Color online) Matrix elements for the sensitivity of 
dose equivalent in water shielding to a given change in each 
partial cross section <Ti y for the 1977 solar minimum GCR 
environment in the thin-shielding limit. 


below the threshold for baryon-antibaryon productions, 
affect the shielded radiation exposure the most [5-7], 
Eq. (9) is a good approximation. 

To study the sensitivity of the shielded dose equivalent 
to fragmentation cross sections, we need first to evaluate 
the change in the dose equivalent caused by the change 
in a partial cross section while keeping all other partial 
cross sections fixed. Therefore, the unitarity constraint 
of Eq. (9) requires us to change the total inelastic cross 
section accordingly, and Eq. (7) becomes 





5<JkjdE.( 10 ) 


For energy-independent changes in the partial cross sec- 
tions, 5<7kj, Eq. (10) reduces to 


5H(x) = px^2 U jk 5a kj , 


u jk =- 


PTPJ 




-Z]Q{Z]L^ + ZlQ{Z*M 


( 11 ) 
L\dE, (12) 


where p represents the density of the shielding material, 
and the relation L k {E) ~ ZlL\(E) has been used for the 
LET spectra with Z k being the charge of particle species 
k. The above matrix element Uj k represents the sensi- 
tivity of the dose equivalent behind a material of areal 
density px to the uncertainty in the partial fragmenta- 
tion cross section a k j- Note that water shielding and 
the 1977 solar minimum GCR environment for elements 
from protons to nickel in the energy range between 0.1 
MeV per nucleon to 1000 GeV per nucleon are used in 
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the calculations unless specified otherwise, and the GCR 
environments in this study are taken from Ref. [8]. 

Fig. 1 shows the matrix elements Ujk as a function 
of projectile and fragment charges. It is evident from 
Fig. 1 that Fe, Si, Mg and O projectiles contribute to 
the highest peaks. For a given projectile, the matrix 
element almost vanishes for either light or heavy (i.e., 
with charges near that of the projectile) fragments while 
peaking for medium-sized fragments. Indeed Eq. (10) 
shows that the contribution from projectile j to 5H(x) 
vanishes at both the light-fragment limit ( Ak —> 0 and 
Zk —■ > 0) and the heavy-fragment limit (k — > j). In a 
simplified case where fragments have the same quality 
factor and the A/ Z ratio as the projectile, we have 

U jk oc -Z}^- +Zt~ Zki-Zj + Z k ), (13) 
A i 


therefore the matrix elements for a given projectile j 
would peak in the middle (at Zk — Zj/'l). We can clearly 
see this feature in the shapes of the fragment distribu- 
tions shown in Fig. 1. 

It is often more useful to evaluate the sensitivity to 
relative changes in the partial cross sections. Considering 
energy-independent relative changes in the partial cross 
sections, Sakj/<Jkj, Eq. (12) reduces to 


5H ( x ) = px S. 


'jk 


j,k 

*= — [jj\-2 
PTpJ L 


fi&kj 

&kj 


,A k 


-Z)Q{Z)U)^+ZlQ{Z^ 


(14) 


\LiakjdE. (15) 



FIG. 2: (Color online) Matrix elements for the sensitivity of 
dose equivalent in water shielding to a given relative change 
in each partial cross section (5crjy/<7jy) for the 1977 solar min- 
imum GCR environment in the thin-shielding limit. 

These results are semi-analytical because they require 
the input of the quality factors and fragmentation cross 


sections. Unlike Ujk of Eq. (12), the shapes of Sjk depend 
on the material due to the Okj term which represents frag- 
mentations in the material. Fig. 2 shows the sensitivity 
matrix elements Sjk with water as the shielding mate- 
rial, where values of the partial cross sections akj(E) are 
taken from Ref. [9]. We see that fragments with Z = 1 
(protons) have the largest effect on the uncertainty of the 
dose equivalent behind shielding, followed by fragments 
with Z = 2 (alpha particles). In the same simplified case 
as for Eq. (13), where fragments have the same quality 
factor and the A/Z ratio as the projectile, we have 

Sjk °c (^Ej + Z’i^j akj — Zki-Zj + Zk)crkj ■ (16) 

Considering Fe projectiles for example, the factor 
Zk{—Zj + Zk) above is the same for fragments with 
Zk = 1 (protons) or fragments with Zk — 25 (Mn), 
but (jpp e (TMnFe- For an Fe projectile at 1.2 GeV 
per nucleon in water shielding, the partial contribu- 
tion to the total inelastic cross section, represented by 
AkGkj/{AjOj ), is about 13.9%, 10.8%, 8.5%, 0.3%, 0.3%, 
1.1%, 8.8%, and 11.9% for neutrons, protons, alpha par- 
ticles, Li, O, Si, Cr and Mn fragments, respectively [9]. 
Thus the cross section to produce a proton fragment, 
crpF e , is about 50 times the value of 0MnFe- Therefore, 
light fragments contribute to the highest peaks in Fig. 2 
partly because of the large partial cross sections to pro- 
duce them. 

Table I lists the top 30 most important cross sections 
for both water and aluminum shielding, ranked accord- 
ing to the absolute values of the corresponding matrix 
elements as given by Eq. (15). We see that the most im- 
portant partial cross sections are for proton and alpha 
particle productions [10] from heavy projectiles such as 
Fe and Si. We also find that in aluminum, another typical 
benchmark material for radiation evaluations, the most 
important cross sections are very similar to those in wa- 
ter shielding. The absolute values of the matrix elements 
in aluminum are smaller in general than those in wa- 
ter shielding, consistent with the fact that aluminum is 
less effective than water as a radiation shielding material 
against GCR. 

Fig. 3 shows one-dimensional reductions of the above 
two-dimensional matrix Sjk to both the projectile and 
the fragment charges. A straightforward summation of 
Sjk over the fragment or the projectile charge would not 
be appropriate because the matrix elements, as well as 
the uncertainties in the partial cross sections, have dif- 
ferent signs in general. In a simplified case where the 
uncertainties in the partial cross sections have the same 
magnitude (\5akj\/&kj = e) but random signs, Eq. (14) 
reduces to 


SH(x) « pxelJ2 Sjk ■ (17) 


Therefore we show 


\Jl2j Sj k as a function of the frag- 
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TABLE I: The top 30 matrix elements in the unit of 
[(cSv/yr)/(g/cm 2 )] ranked according to their absolute values 
for the 1977 solar minimum GCR environment. 



For water shielding 


For aluminum shielding 

j 

k 

Sjk 

j 

k 

—Sjk 

26 

1 

0.2800 

26 

1 

0.2204 

26 

2 

0.2488 

26 

2 

0.1691 

14 

1 

0.1639 

14 

1 

0.1173 

8 

1 

0.1495 

8 

1 

0.1019 

12 

1 

0.1280 

12 

1 

0.0883 

14 

2 

0.0982 

14 

2 

0.0634 

12 

2 

0.0647 

16 

1 

0.0395 

8 

2 

0.0627 

12 

2 

0.0388 

10 

1 

0.0545 

10 

1 

0.0380 

16 

1 

0.0507 

20 

1 

0.0375 

6 

1 

0.0490 

6 

1 

0.0348 

20 

1 

0.0483 

22 

1 

0.0319 

8 

7 

0.0447 

24 

1 

0.0315 

22 

1 

0.0422 

8 

2 

0.0277 

8 

6 

0.0420 

25 

1 

0.0239 

24 

1 

0.0405 

24 

2 

0.0235 

14 

12 

0.0373 

22 

2 

0.0226 

16 

2 

0.0348 

20 

2 

0.0201 

24 

2 

0.0345 

8 

7 

0.0200 

20 

2 

0.0342 

13 

% 

0.0199 

14 

10 

0.0331 

18 

i 

0.0193 

22 

2 

0.0325 

25 

2 

0.0185 

14 

11 

0.0320 

16 

2 

0.0180 

8 

5 

0.0314 

8 

6 

0.0172 

25 

1 

0.0306 

8 

5 

0.0166 

12 

11 

0.0296 

19 

1 

0.0156 

13 

1 

0.0287 

23 

1 

0.0155 

14 

9 

0.0272 

7 

1 

0.0155 

12 

10 

0.0271 

14 

11 

0.0146 

25 

2 

0.0269 

14 

10 

0.0136 


ment charge Z k in Fig. 3(a) and ^jYlk $ jk 88 a function 
of the projectile charge Zj in Fig. 3(b). Fig. 3(a) clearly 
shows the dominance of fragments with Z = 1 (protons) 
and Z = 2 (alpha particles) while Fig. 3(b) shows that 
fragmentations of projectiles such as Fe, Si, Mg and O 
have the largest effects on the uncertainty of the dose 
equivalent behind shielding. 


III. BEYOND THE THIN-SHIELDING LIMIT 


Eqs. (14) and (15) derived in the thin-shielding limit 
can be written as 


clSH(x) 

dx 


p 'y , Sjk (p 

:i,k 


J fi&kj 
<Jkj 


(18) 



FIG. 3: Fragment (a) and projectile (b) distributions for the 
sensitivity of dose equivalent to a given relative change in the 
partial cross sections. 
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FIG. 4: (Color online) Matrix elements for the sensitivity 
of dose equivalent at 20 g/cm 2 in water shielding to a given 
relative change in each partial cross section. 


S jk (x) = — / J,[E,x 
PTP J 

+ ZlQ(ZlLr)} LiakjdE 


Z jQ(ZjLi)~ 


(19) 


where Sjk{x) represents the sensitivity matrix element 
at any given depth x in a shielding material. Note that 
Eq. (18) gives the rate of the dose equivalent change, i.e., 
the dose equivalent change per. unit depth, at depth x 
from energy-independent relative changes in the partial 
cross sections. 
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FIG. 5: Fragment (a) and projectile (b) distributions for the 
sensitivity of dose equivalent at 20 g/cm 2 in water shielding 
to a given relative change in the partial cross sections. 


As an example, using the deterministic radiation trans- 
port code HZETRN [8] we have calculated the particle 
spectra at 20 g/cm 2 in water shielding, which were then 
used as the input Jj to evaluate the sensitivity matrix 
elements according to Eq. (19). Fig. 4 shows the sen- 
sitivity matrix Sjk at 20 g/cm 2 of water shielding, and 
Fig. 5 shows the one-dimensional reductions of that sen- 
sitivity matrix. First of all, compared with the results 
in the thin-shielding limit (i.e, at x ~ 0 g/cm 2 of water 
shielding), we see that at 20 g/cm 2 of water shielding pro- 
jectiles such as O and Mg are relatively more important 
while the Fe projectile is relatively less important. This 
is mainly because heavier projectiles have larger inelastic 
cross sections, and thus their fluences tend to decrease 
more inside materials. Secondly, Fig. 5(a) shows that 
the fragment distribution at this depth has almost the 
same shape as that in Fig. 3(a), therefore cross sections of 
light fragment productions are still the dominant source 
of uncertainty. Also, we see that the sensitivity matrix 
elements at 20 g/cm 2 of water shielding are much smaller 
in magnitudes, which reAects the overall decrease of the 
dose equivalent behind shielding as well as the decreased 
efficiency of fragmentations after a significant amount of 
shielding. 

The total sensitivity of the dose equivalent behind a 
shielding material of thickness L to cross sections can be 
evaluated according to Eq. (18) as 



dSH(x) , 

— ax 

ax 



Sjk{x)dx 


6(7 fcj 
&kj 


( 20 ) 


In a simplified case where the errors in the partial cross 
sections have the same magnitude as well as the same 
sign ( 5a kj /a kj = e'), Eq. (20) gives pe' £ jVc f Q L S jk {x)dx. 


Conversely, in the same simplified case as for Eq. (17), 
where the uncertainties in the partial cross sections have 
the same magnitude (Idcpyl/cqy = e) but random signs, 
Eq. (20) reduces to 



dSHix) , 
— — ax 

ax 



(21) 


In any case, the sensitivity for a thick shielding is dom- 
inated by the depths at which the sensitivity matrix el- 
ements have the largest values. Fig. 6 show s the projec- 
tile distributions, represented by \jYlk ‘S'jfo at f° ur dif- 
ferent depths in water shielding, and it shows that the 
sensitivity matrix elements for charged fragments have 
larger values at smaller depths for GCR environments 
[10], This suggests that the semi-analytical results in the 
thin-shielding limit are representative of the total sen- 
sitivity for GCR propagation through a relatively thick 
material. 



Projectile Z ( 


FIG. 6: (Color online) Projectile distributions for the sensi- 
tivity of dose equivalent at 0, 5, 10 and 20 g/cm 2 in water 
shielding to a given relative change in the partial cross sec- 
tions. 


IV. DISCUSSIONS 

To illustrate the effects of the unitarity constraint of 
Eq. (9), we show in Fig. 7 the matrix elements without 
imposing the unitarity relation. In this case the terms 
with a negative sign in Eqs. (7), (10), (12), (15) and 
(19) have been neglected, leading to the following matrix 
elements: 

Sf k = — [ df/lQ{ZlL x )L,o k ,dE. ( 22 ) 

PTP J 

In Fig. 7 we find that heavy fragments with charges near 
that of the projectile are the most important for a given 
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projectile; an expected feature due to the strong Zk de- 
pendence in Sjj.. On the other hand, results with the 
unitarity constraint shown in Fig. 2 show totally different 
fragment distributions. We also see that the matrix ele- 
ments without the unitarity constraint in Fig. 7 are pos- 
itive, while those in Fig. 2 with the unitarity constraint 
are mostly negative, indicating that in most cases an in- 
dependent increase of a partial cross section decreases 
the dose equivalent behind shielding. 
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FIG. 7: (Color online) Matrix elements, without imposing the 
unitarity constraint, for the sensitivity of dose equivalent to 
a given relative change in each partial cross section. 

The semi-analytical method provides the theoretical 
understanding of the sensitivity of the shielded dose 
equivalent to uncertainties in the partial cross sections. 
The next step is to apply the semi-analytical method in 
combination with our knowledge on fragmentation cross 
sections and on their uncertainties, where other imple- 
mentations of the unitarity constraint may be needed. 
After a particular Okj is changed, instead of changing 
the total inelastic cross section crj to satisfy Eq. (9) a 
combination of the total inelastic cross section and one 
or more other partial cross sections may be changed. In 
this case, however, there is no unique definition for the 
sensitivity to a single partial cross section because the 
effect will always correspond to changes in multiple par- 
tial cross sections [11]. As our goal is to establish the 
priorities for partial cross sections to be further studied, 
how to implement the unitarity constraint should depend 
on our knowledge of these cross sections, especially the 
sizes of their errors and the correlations among them. 
Correlations in experimental data are certainly possible 
among partial cross sections, as well as between partial 
cross sections and the total inelastic cross section. 

To explore other ways of implementing the unitarity 
constraint, we consider two examples where total inelas- 
tic cross sections crj are kept the same. The first example 




FIG. 8: (Color online) Correlated matrix elements for the 
sensitivity of dose equivalent to a given relative change in 
each partial cross section a kj, with the unitarity constraint 
of Eq. (9) being satisfied by fixing oy- and changing all other 
partial cross sections by the same fraction. 



FIG. 9: (Color online) Correlated matrix elements for the 
sensitivity of dose equivalent to a given relative change in 
each partial cross section cqy , with the unitarity constraint of 
Eq. (9) being satisfied by fixing cxj and changing 


assumes that the total inelastic cross sections are known 
exactly but uncertainties in the partial cross sections are 
mostly uncorrelated. Thus a change in Ukj can be ac- 
companied by changes in all other partial cross sections, 
(7y(i ^ k), by the same fraction in order to satisfy the 
unitarity constraint of Eq. (9). The resulting sensitivity 
matrix elements are shown in Fig. 8, and we see that 
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both heavy and light fragments are important in this 
case. The second example also assumes that the total 
inelastic cross sections are known exactly but there are 
misidentifications of particles with charge Zk as particles 
with charge Zk - 1 (for all k), where an under-measured 
< Tkj data would be correlated with an correspondingly 
over-measured <rk-i,j data. In this case a change in 
is accompanied by a corresponding change in crk-i,j in 
order to satisfy Eq. (9). The resulting sensitivity ma- 
trix elements are shown in Fig. 9, and we see that frag- 
ments with charges near that of the projectile for certain 
medium-sized projectiles are the most important. This 
feature can be understood using the matrix elements: 


pTpJ L Ak-x. 


k/iakjdE,(23) 


which, in the simplified case where fragments have the 
same quality factor and the A/Z ratio as those of the 
projectile, result in S jjj. oc SK/Zk- 

We emphasize that a smaller relative uncertainty in the 
total inelastic cross section than in the partial cross sec- 
tions does not necessarily mean that one has to keep the 
total inelastic cross section aj fixed, because the smaller 
uncertainty may simply result from the fact that the to- 
tal inelastic cross section is the sum of the partial cross 
sections [5, 7]. Consider the simplified case where there 
are Nj partial cross sections for projectile j with similar 
values of Ak<?kj and similar magnitudes of uncertainty 
(\5okj\/<rkj — e) with random signs, while the other par- 
tial cross sections have negligible values of Akcrkj, Eq.(9) 
then gives 


5(Xj e 


(24) 


therefore the uncertainty in the total inelastic cross sec- 
tion is much smaller. 

To see the extent to which the sensitivity matrix de- 
pends on solar cycle variations of the GCR environments, 
Fig. 10 shows the one-dimensional reductions of the sen- 
sitivity matrix Sjk in water shielding for the 1989 solar 
maximum GCR environment. Compared with results for 
the 1977 solar minimum GCR environment as shown in 
Fig. 3, we find that the sensitivity matrix has similar 
features but the overall normalization is lower mainly 
because of the smaller fluences of GCR particles during 
solar maximum. 

This study addresses the important cross sections for 
the protection against the exposure to galactic cosmic 
rays. Such a study on the protection against solar parti- 
cle events will be useful since they also pose a serious radi- 
ation hazard. However, we note that the straight-ahead 
approximation assumed in this semi-analytical method 
does not apply well to SEPs because of their dominant 
low-energy part of the spectra, and charge conservation 
may impose an important constraint as the majority of 
the nuclear fragmentation processes will occur below the 
pion production threshold. 



FIG. 10: Fragment (a) and projectile (b) distributions for 
the sensitivity of dose equivalent to a given relative change in 
the partial cross sections for the 1989 solar maximum GCR 
environment. 


V. CONCLUSION 

In this study we have developed a semi-analytical 
method to calculate the contribution to the uncertainties 
of the dose equivalent after shielding from each partial 
cross sections of nuclear fragmentation in the shielding 
materials. This method is able to point to the most im- 
portant partial cross sections of nuclear fragmentations 
for radiation protection of astronauts against galactic 
cosmic rays in long-term space missions. By varying each 
partial cross section independently while satisfying the 
unitarity relation due to the baryon number conserva- 
tion, we find that light fragment productions from heavy 
ions such as Fe, Si, Mg and O appear the most important. 
Thus, better theoretical and experimental understanding 
of these cross sections will help reduce the uncertainty in 
the radiation health risk predictions in long-term space 
missions. The unitarity relation required by baryon num- 
ber conservation is shown to have significant effects on 
the sensitivity matrix as it correlates the uncertainties 
in the cross sections. At a finite depth in the shielding 
material, projectiles such as O and Mg become relatively 
more important while the Fe projectile becomes relatively 
less important, and the matrix elements have much lower 
magnitudes. Therefore, the thin-shielding results should 
be representative of the total sensitivity for thick shield- 
ing. 

The next step is to apply the semi-analytical method in 
combination with our knowledge on fragmentation cross 
sections and their uncertainties. We have provided two 
simple examples, where different implementations of the 
unitarity relation may be used depending on current un- 
certainties in our knowledge of the fragmentation cross 
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sections. We have also shown that the sensitivity matrix 
for a solar maximum GCR environment is similar to that 
for a solar minimum GCR environment except that the 
overall magnitude is much smaller. 
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